Free energy cascade in gyrokinetic turbulence 
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In gyrokinetic theory, the quadratic nonlinearity is known to play an important role in the dynam- 
ics by redistributing (in a conservative fashion) the free energy between the various active scales. In 
the present study, the free energy transfer is analyzed for the case of ion temperature gradient driven 
turbulence. It is shown that it shares many properties with the energy transfer in fluid turbulence. 
In particular, one finds a forward (from large to small scales), extremely local, and self-similar cas- 
cade of free energy in the plane perpendicular to the background magnetic field. These findings 
shed light on some fundamental properties of plasma turbulence, and encourage the development of 
large eddy simulation techniques for gyrokinetics. 



Fully developed turbulence is fundamentally linked to 
a conservative transfer of (free) energy in wavenumber 
space from drive to dissipation scales [l[. While the re- 
spective cascade dynamics for simple fluids (described 
by the Navier-Stokes equation) has been the subject of 
countless studies and is fairly well understood, the situa- 
tion is quite different for turbulent plasmas, both at large 
scales (compared to the gyroradii of the particles) - de- 
scribed in the context of magnctohydrodynamics - and, 
in particular, at small scales - described by the gyroki- 
netic equations Q. The latter case, in which one deals 
with a gyrocenter distribution function in three spatial 
dimensions as well as two velocity space dimensions (the 
third velocity space coordinate can be removed analyti- 
cally in a low- frequency ordering), shall be the focus of 
the present work. 

In three-dimensional Navier-Stokes turbulence, the ki- 
netic energy is conserved by the convective nonlinearity. 
It is usually assumed to be injected into the system at the 
largest scales through mechanical forcing, and to be dis- 
sipated at the smallest scales by viscous effects. The role 
of the nonlinearity is then to transfer the kinetic energy 
from the large scales to the small ones in what is usu- 
ally referred to as a cascade process. In the gyrokinetic 
formalism, on the other hand, the free energy acts as the 
quadratic conserved quantity (see, e.g., Ref. Q and var- 
ious references therein). It is usually injected into the 
system at large scales via the background density and 
temperature gradients, and expected to be dissipated at 
small (space and/or velocity space) scales. It is antici- 
pated that one role of the nonlinear term in gyrokinetic 
turbulence is to transfer the free energy from the largest 
perpendicular scales to the smallest ones 0-0, but a 
definitive investigation of the free energy transfer dynam- 
ics in a self-driven, three-dimensional system (which is 
the standard case for magnetically confined plasmas) is 
still lacking and shall be provided for the first time in the 
present Letter. 



Our study is based on numerical solutions of the non- 
linear gyrokinetic equations obtained by means of the 
Gene code 0-0- Although Gene is able to treat an 
arbitrary number of fully gyrokinetic particle species as 
well as general toroidal geometry, magnetic field fluctu- 
ations, and collisions, these features shall not be used 
here. Instead, we will focus on the reduced problem of a 
single ion species, adiabatic electrons, electrostatic fluc- 
tuations, and a large aspect-ratio, circular cross-section 
model equilibrium. For this simplified case, the respec- 
tive (appropriately normalized) equations read (for de- 
tails, see Ref. @): 
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Here, the total distribution function Fj of species j is 
split into a Maxwellian part F oj = 7r _3 / 2 e _( V+' lBo ) 
and a perturbed part fj, and the nonadiabatic part 
of fj is given by hj = fj + (qj<j>i/T j) F j where <j> x 
is the gyro-averaged electrostatic potential, hj and fj 
depend on the gyrocenter position r = (x,y,z), the 
parallel velocity vu , the magnetic moment /i, and the 
time t. As indicated already above, all simulations in 
this paper are performed in s — a geometry [Io| with 
a = 0, for which the curvature terms are given by 
JC X = — 2sinz and IC y = — 2 (cos 2; + sz sin z). Fur- 
thermore, VTj = {ZTjo/rrij) 1 / 2 is the thermal velocity, 
u) n . = —Rd\ogriQj/dx and wt, = —RdlogToj/dx are 
the normalized background density and temperature gra- 
dients, nij and qj are the mass and charge of species j. 
The equilibrium magnetic field is taken to be B = Bq B vc { 
where B re { is the reference magnetic field on the magnetic 
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axis. Finally, the Poisson brackets are defined by 
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Note that in Eq. ([1}, the second term is responsible for 
the injection of free energy into the system. The third 
through fifth terms are, respectively, the curvature, non- 
linear, and parallel terms, none of which acts as a source 
or sink of free energy. Since the simulations presented 
below are done without collision operator, the numerical 
scheme used in Gene is not dissipative, and a statisti- 
cal steady state cannot be reached without some form of 
dissipation Here, hyperdiffusion terms T> z and V v ^ 
are added to dissipate fine-scale fluctuations in z and v\\ 
(for details, see Ref. flij). 

Eq. ((T|) is complemented by the gyrokinetic Poisson 
equation which is used to determine the self-consistent 
electrostatic potential: 
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Here, Jo is the Bessel function and To(bj) = e~ bj Io(bj) 
with the modified Bessel function Iq. The (dimension- 
less) arguments bj and Xj are defined, respectively, as 



and the (positive definite) dissipative term 

V = -Hj dA l^ h i (P*h + V n fi) ■ (8) 
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The quantity £ plays the same role in gyrokinetic turbu- 
lence as the kinetic energy in fluid turbulence Q. 

The transfer of free energy between different modes in 
the saturated turbulent state is induced by the nonlin- 
ear term. Although it does not affect the global value 
of the free energy (numerically, this is satisfied in Gene 
up to machine precision), it can change, e.g., the value 
of this quantity associated with particular perpendicular 
wavenumbers. Following the procedure used for study- 
ing energy transfer in Navier-Stokes and in MHD turbu- 
lence Ha-ai, 

we decompose the perpendicular wavevec- 
tor plane into domains and measure the free energy trans- 
fer between these domains. The set of domains {de} is 
assumed to be a partition (no intersection between the 
domains and all domains together cover the entire plane) . 
The distribution function and electrostatic potential can 
then be written as a sum over all contributions for which 
the perpendicular wavevectors lie in the domain dg. As 
a consequence of the Parseval theorem, the free energy 
can also be split into parts which are associated to the 
domains dp: 
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where Oj = (qjB a )/ (rrijc) and k± is the perpendicular 
wave number. 

In the absence of drive and dissipation, the gyrokinetic 
equations, Eqs. ([1} and are known to conserve the 
free energy £ (see, e.g., Refs. @,Q) which is usually split 
into two quadratic parts according to £ = £f + £$ with 
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Here, J dA = j d 3 x j 7r_B n 0i dv\\dfi denotes phase-space 
integration. The evolution equation for the free energy 
is given by 
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In the problem considered hereafter, both the entropy 
and electrostatic contributions to the free energy are con- 
served separately by the nonlinearity Af. It is thus legiti- 
mate to consider the entropy conservation independently 
from the conservation of the electrostatic energy. The 
evolution of £j due to the nonlinear term can be ex- 
pressed as 
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where we have used the property J dAfjfj = J dAfjfj 
which is easily proven and expresses the fact that the 
contributions fj are orthogonal "vectors" if their scalar 
product is defined as the integration over A. Introducing 
the explicit form of the nonlinearity, one obtains 
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in terms of the source term 



where the three-domain interaction terms are defined as 
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ff\ Eq. (fTTj) shows that the evolution of the entropy associ- 
ated to the domain di is the sum of triple interactions 



3 



between wave vectors associated to the domains dg, dg 1 
and dg 2 . This is not a surprise since, like in the Navier- 
Stokes equation, the quadratic nonlinearity in the gy- 
rokinetic equation is responsible for triadic interactions 
between the Fourier modes. Proposing a clean defini- 
tion of the energy transfer between two domains might 
thus be problematic in such a picture. However, con- 
sidering the structure of these three-domain interaction 
terms, the following two-domain interaction terms is a 
natural quantity to investigate: 
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These two-domain interaction terms will be interpreted 
as the energy transfers between the domains dg and dgi, 
even if the redistribution of the free energy between the 
different domains by the nonlinear term cannot be fully 
understood without considering triadic interactions. As 
a consequence of the Poisson bracket structure, it is easy 
to show that T^' e = —T* ' e , which reinforces the inter- 
pretation in terms of free energy exchange. Indeed, if the 
domain dg is considered to receive a certain amount of 
free energy per unit of time Ty from the domain dgi, 
then the domain dgi is seen as loosing exactly the same 
amount of free energy per unit of time in profit of the 
domain dg. The same approach can be used to define 
three-domain and two-domain interaction terms for the 
electrostatic part of the free energy with the following 
definitions: 
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The complete dynamical equation for £ l then reads 

¥ = E T f +E T f +0*-* ( 16 ) 

m v g> 

where the source and dissipation terms, Q e and T> 1 
given, respectively, by Eqs. ^ and ©, using h~j, /|. 



are 
and 



The free energy transfer terms defined above are now 
evaluated from a numerical simulation using Gene. The 
physical parameters employed in this context correspond 
to a widely used standard case of collisionless ion temper- 
ature gradient (ITG) turbulence known as the Cyclone 
Base Case [l7|. The simulation domain is about 125 
ion gyroradii wide in the perpendicular directions, and 
256 x 64 x 64 x 32 x 8 grid points are used in (x, y, z, vu , /i) 
space. For further analysis, the perpendicular wavevec- 
tor plane is divided into shells dg = {kj_ such as Kg < 
kjj < Kg + i} where the shell boundaries Kg are cho- 
sen to grow algebraically Kg+i = \Kg, with A = 2 1 /' 5 



between shell £ = 3 and £ = 24. The first shell bound- 
aries have been chosen differently (K\ = 0, if 2 = 0.2, 
K 3 = 0.3) in order to ensure that enough modes belong 
to those shells. Moreover, in order to limit the number 
of shells, the last shell (£ = 25) is wider and limited by 
K25 = 6.3 and K 26 = |kj max = 14.6. 
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FIG. 1. Shell decompositions in perpendicular wavenumber 
space of the drive (Q e ) and dissipation (T> ( ) terms (as well as 
their sum) from a Gene simulation of ITG turbulence. 

Fig. [T] shows the numerical results for the source and 
dissipation terms (averaged over time during the satu- 
rated phase of the simulation). As expected, the injec- 
tion of free energy is well localized at low k± . However, 
as it turns out, the dissipative terms are not just ac- 
tive in the high k±_ range, but throughout the entire kj_ 
spectrum, including the drive range. An explanation of 
this phenomenon may be provided in terms of the non- 
linear coupling to damped eigenmodes, as is discussed in 
Ref . [l8| ■ There is a net source of free energy up to shell 
I = 9 and a net dissipation beyond that. The peak in the 
dissipation may be due to the fact that the largest shells 
are not complete because we are using a discretization in 
and k y . 

The corresponding shell-to-shell free energy transfer 
terms are shown in Fig [31 and various interesting fea- 
tures can be observed there. First, the entropy transfer 
is larger than the electrostatic energy transfer by almost 
two orders of magnitude, thus dominating the total free 
energy transfer. Second, the free energy transfer is from 
the large scales to the small ones (as one might have ex- 
pected); the transfer is systematically negative for I' > I 
and, due to the antisymmetry property, systematically 
positive otherwise. Third, the free energy transfer is very 
local in wavenumber space. Indeed, only values of T t ^ t 
with £ close to £' are significantly different from zero. In 
practice, for \£ — 1'\ > 5 the free energy transfers almost 
vanish. This corresponds to a ratio of wave numbers be- 
tween the two shells of the order of two. Fourth, a lim- 
ited self-similarity range can be identified for £ between 
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FIG. 2. Shell-to-shell transfer in perpendicular wavenumber 
space of entropy (a) and electrostatic energy (b) from a Gene 
simulation of ITG turbulence. 

13 and 20. Indeed, in this range, the total transfers T t ^ t 
seem to depend on I — £' only, and not on the two in- 
dices separately. Considering the limited resolution of 
the simulation analysed here, this property is rather un- 
expected. Indeed, the analysis of the source and dissipa- 
tion terms (see Fig. [J) does not show the existence of a 
range of scales in which both these terms would be neg- 
ligible. However, they are obviously sufficiently small to 
allow for cascade dynamics to develop (see also Ref . [l8| ) ■ 
This is actually to be expected if the nonlinear frequen- 
cies characterizing the free energy transfer exceed the 
linear ones characterizing the dissipation. 

Interestingly, the spectral transfer of free energy in 
gyrokinetic turbulence thus exhibits various similarities 
with respect to the kinetic energy transfer measured in 
fully developed Navier-Stokes turbulence, although this 
is not all clear a priori. In particular, there is a (strongly) 
local, self-similar forward cascade - despite the absence 
of an inertial range. Insights like these may be expected 
to guide the application of large-eddy simulation tech- 



niques [2(| to gyrokinetics. Here, the idea is to only 
retain the dynamics of the largest scales while the small- 
est ones are modelled. Indeed, if the smallest scales are 
proven to act systematically as a sink of free energy like 
it was the case here, it is reasonable to propose a dissipa- 
tive model for these small scales and consequently reduce 
as much as possible the numerical resolution. On such 
a basis, it may well become possible to reduce the com- 
putational effort for gyrokinetic turbulence simulations 
by a significant amount. The present work represents a 
relevant step in that direction. 
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